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ABSTRACT 

The magnetic Prandtl number PrM is the ratio of viscosity to resistivity. In 
astrophysical disks the diffusion of angular momentum (viscosity) and magnetic 
fields (resistivity) are controlled by turbulence. Phenomenological models of 
the evolution of large scale poloidal magnetic fields in disks suggest that the 
turbulent magnetic Prandtl number PrM,T controls the rate of escape of vertical 
field from the disk; for PrM,T ^ R/H vertical field diffuses outward before it can 
be advected inward by accretion. Here we measure field diffusion and angular 
momentum transport due to MHD turbulence in a shearing box, and thus PrM,T, 
by studying the evolution of a sinusoidal perturbation in the magnetic field that 
is injected into a turbulent background. We show that the perturbation is always 
stable, decays approximately exponentially, has decay rate cx fc^, and that the 
implied PrM,T ~ 1- 

Subject headings: accretion, accretion disks, magnetohydrodynamics 



1. Introduction 

Astrophysical disk evolution may be controlled in part by magnetic fields that are coher- 
ent over scales of order the radius R. Large scale fields are an e ssential element of theoret ical 



models for the laun ching and coUimatio n of disk winds (e.g. iBlandford fc Paynelll982l : for 



a recent review see 



Pudritz et al.l 120071 ). which lead to disk evolution because they exert 
direct torques on the surface of the disk. Large scale fields may also control the strength 
of turbulent angular momentum diffusion (in dimensionless form, a) within the disk, since 
numerical experiments in unstratified "shearing boxes" suggest that a is proportional to the 
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mean field strength (e.g. iHawley et al.lll995l ). So what determines the strength of the large 



scale field in well ionized disks? 

A few of the physical processes affecting the large scale field strength in thin disks can be 
easily listed: (1) advection of flux by large scale flows in the disk (e.g., inward advection by 
accretion, but also possibly meridional circulation); (2) turbulent diffusion; (3) dynamical 
interaction between large scale fields and large scale flows (e.g. enhanced accretion due 
to external torques); (4) control of small scale MHD turbulence by large scale fields (e.g. 
through the influence of a mean field on a); (5) generation of large scale fields by small scale 
MHD turbulence (a large scale dynamo); (6) introduction or removal of magnetic flux at 
the disk boundaries. These processes are difficult to study numerically because they involve 
nonsteady flows and a large dynamic range in length scale [R : scale height H) and time 
scale (viscous timescale : dynamical timescale). 

A starting point for understanding large scale field evolution is the phenomenological 
model of Van Ballegooijen (1989; hereafter VB89). He considers a passive large scale field 
advected inward by accretion and diffused by turbulence in the disk (processes [1] and [2] 
above). Turbulence is modeled using a turbulent viscosity i^t and turbulent resistivity ?7t. 
The evolution of the poloidal field is then governed by the induction equation in vector 
potential form: 

dtA^ = mdRl^dniRA^)] + T^^dlA^ + VR^dniRA^), (1) 

ri K 

where A^ is the azimuthal component of the vector potential. A^ labels field lines; when 
dtA^ 7^ the field lines move radially through the disk with speed ~ dtA(f,/dRA^. In the 
limit that rj^ a. vertical field would be advected inward by accretion and vertical field 
strength would increase with time. In the limit that z/ the field lines simply diffuse out 
of the disk. Where do advection and diffusion balance? 

VB89 give a surprising answer that can be understood as follows, assuming the field 
and disk are symmetric about the midplane. At the midplane the first term in ([T]) (radial 
diffusion of vertical field) is ~ titBz/R, where = {1/ R)dR{RA^) and we assume that 
Or ^ 1/ R. The second term (vertical diffusion of radial field, which can nevertheless cause 
radial motion of field lines) depends on the field geometry above and below the disk. If we 
as sume the field lines enter a nd exit the disk at an angle of order unity (as in the wind model 



of iBlandford fc Payndll982l ). Br = —d^A^ ~ =f_B^ at z = ±H and Br = at the midplane, 
by symmetry. Then rj^dlA^ ~ rj^Bz/H. Using the usual viscous disk estimate vr ~ i^/R, 
the final term (field advection) is of order vB^/R. To summarize, the terms on the right 
hand side of ([1]) are in the ratio rj-j^Bz/R : rj^Bz/H : vBz/R. Evidently the first term is 
negligible in comparison to the second for a thin disk, provided the turbulent diffusion can 
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be described by a scalar diffusion coefficient. The second and third third can balance when 
i^t/ Vt ^ R/H] diffusion and advection balance when PrM,T ~ R/H. 



It is plausible that in disk turbulence PrM,T ~ 1 (e.g. lYousef et al.l (120031 )). Then 
outward field diffusion occurs on a timescale RH/ v. Large scale poloidal fields would vanish 
from the disk absent a dynamo that regenerates the field on the s ame timescale (process 
[5] above). Similar conclusions have been reached by iLubow et al.l (119941) in the cori t ext o f 
magnetically generated outflow/jet models. More complex models by lHeyvaerts et al.l (119961 ) 
also support such a picture. 

There are of course ways to avoid the loss of large-scale field implied by the VB89 
mode l, which is based on a pu rely phenomenological model for evolution of the disk and 
field. ISpruit &: UzdenskyI (120051 ) discuss a model that reduces turbulent diffusion by grouping 
large scale vertical magnetic fields into bundles in the disk through flux expulsion. In these 
bundles the fields are strong enough to quench turbulence an d thus avoid outward d if fusion 
(processes [3] and [4] above). It has also been suggested by lUzdensky fc GoodmanI (120081 ) 
that disk atmospheres might develop loop-like large scale coro nal structures that delocalize 
disk evolution by transmitting angular momentum and energy. iRothstein fc Lovelacd (120081 ) 
consider the possibility that field diffusion is suppressed by a coronal layer where f a > c^, so 
that the MRI is suppressed and t^t is reduced. Another possibility is that disk evolution is 
driven by external torques associated with an MHD (Blandford-Payne type) wind running 
along the poloidal field lines. If the field is strong enough then the inflow speed vr may be 
large enough to compete with outward diffusion of field lines even when PrM,T ~ 1 (process 
[3] above). [J A final possibility is that PrM,T ^ 1- 

In this paper we measure the turbulent magnetic Prandtl number Prn.T directly from 
shearing box simulations. We infer the turbulent viscosity z/t from the turbulent shear 
stress Wxy^T (in dimensionless form, a), which controls diffusive radial transport of angular 
momentum in disks. We infer the turbulent resistivity t^t by tracking the evolution of a 
sinusoidal disturbance in the magnetic field that is imposed on an already turbulent state. 

It is worth emphasizing that the turbulent Prandtl number PrM,T is funda mentally dif- 



ferent from the Prandtl number Prjv/ associated with microscopic processes (e.g.. lBalbus fc Henri 



^Although a naive estimate suggests this will not work. Suppose that the external torque per unit 
area r ~ RM^^ ^ {Y./ H)R{va)'^ , where Mij = magnetic stress tensor, (va) = the midplane Alfven speed 
associated with the large scale field and S = surface density. Then the induced Vr ^ T/{flJ^R) ^ {va)'^/cs, 
so in eq.(l) the second and third terms balance for 7/t ^ (va)"^/^- If ?7t ^ {5v\)/n {5va = the Alfven 
speed associated with disk turbulence), then (diffusion/advection) ~ {5v\) / {va)^ ■ Shearing box experiments 
suggest that this is large compared to 1 unless the poloidal field is so strong that turbulence is suppressed. 
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2008 ). Recent numerical experiments (IFromang &: Papaloizoul 120071 : iLesur fc Longaretti 



20071 ) suggest that Ptm can influence the saturation level of MRI driven disk turbulence 



at low Reynolds number. 

The paper is organized as follows. In §2 we give a simple description of the local model 
and summarize our numerical algorithm. In §3 we describe the numerical procedure for 
measuring ?7t. We report ?7t and discuss its dependence on the model parameters. §4 We 
explain how we calculate PrM,T and discuss our results. 



2. Local Model and Numerical Methods 



Our starting point is the local model for disks. It is obtained by expanding the equations 
of motion around a circular-orbiting coordinate origin at cylindrical coordinates (r, 0, z) = 
(ro, Qot + (po, 0), assuming that the peculiar velocities are comparable to the sound speed and 
that the sound speed is small compared to the orbital speed. The local Cartesian coordinates 
are obtained from cylindrical coordinates via {x, y, z) = {r—Vo, ro[(f)—flot—(f>o], z). We assume 
throughout that the disk is isothermal (p = cfp, where Cg is constant), and that the disk 
orbits in a Keplerian (1/r) potential. 



In the local model the momentum equation of ideal MHD becomes 



dv „ 

— + V ■ Vv + c: 

ot 



2nxv- 3n'^xx = 0. 



(2) 



p Svrp 47rp 

The final two terms in equation ([2]) represent the Coriolis and tidal forces in the local frame. 
Our model is unstratified, which means that the vertical gravitational acceleration —Q'^z 
usually present in Keplerian disks is ignored. The box has size L^^ x Ly x L^. 

Our model c ontains no explicit dissipation coefficients. Recent mo dels with explicit 
scalar dissipation (IFromang fc Papaloizoul 120071 . iLesur &: Longarettil 120071 ) have shown that 



the saturated field strength in magnetized disk turbulence depends on the viscosity u and 
resistivity rj, and that ZEUS has an effective magnetic Prandtl number PrM,T = '^/v ^ ^ 
(IFromang fc Papaloizoul 120071 ). 



The orbital velocity in the local model is 



^orb 



-Qxij. 
2 ^ 



(3) 



This velocity, along with a constant density and zero magnetic field, is a steady-state solution 
to equation ([2]). If the computational domain extends to |x| > {2/3)H = {2/3)cs/Q, then 
the orbital speed is supersonic with respect to the grid. 
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The loca l model can be stud ied numerically using the "shearing box" boundary con- 
ditions (e.g. iHawley et al.l Il995l ). The boundary conditions on the y boundaries of the 
box are periodic, while the x boundaries are "nearly periodic", i.e. they connect the radial 
boundaries in a time-dependent way that enforces the mean shear flow. We also use periodic 
boundary conditions in the vertical direction; this is the simplest possible version of the 
shearing box model. 



Our models are evolved using ZEUS (IStone fc Normanlll992l ). ZEUS is an operator-split, 
finite difference scheme on a staggered mesh. It uses artificial viscosity to capture shocks. 
B For the magnetic field evolution ZEUS uses the Method of Characteristics-Constrained 
Transport (MOC-CT) scheme, which is designed to accurately evolve Alfven waves (MOC) 
and to preserve the V • -B = constraint to machine precision (CT). 

We have modif i ed ZE US to include "orbital advection" (Massetl 2000 : Gammie 2001 



Johnson fc Gammid l2005l ) with a magnetic field (j Johnson et al.ll2008l ). Advection by the 



orbital component of the velocity Vorb (which may be supersonic with respect to the grid) is 
done using interpolation. With this modification the timestep condition At < CAx/ {\5v\ + 
Cmax) {cmax = maximum wave speed and C = Courant number) depends only on the per- 
turbed velocity 6v = v — Vorb rather than v. So when \vorb\ ^ Cmax (for shearing box 
models with f^/c^ < 1, when L > H) the timestep can be larger with orbital advection, and 
computational efficiency is improved. 

Orbital advection also improves accuracy. ZEUS, like most Eulerian schemes, has a trun- 
cation error that increases as the speed of the fluid increases in the grid frame. In the shearing 
box without orbital advection the truncation error would then increase monotonically with 

Orbital advection reduces the arnplitu de of the truncation error and also makes it more 
nearly uniform in (j Johnson et al.ll2008l ). 



We have also implemented an additional procedure to remove the radi ally dependent 



numer ical dissipation in large shearing box simulations that was reported in iJohnson et al. 
(120081 ). We do so by systematically shifting the entire box by a few grid points in the 



radial direction at t = n2Ly/ (3f2L^, n = 1, 2, 3, . . . (this is when the shearing box boundary 
conditions are exactly periodic). After the shift we execute a divergence cleaning procedure 
to remove the monopoles that build up due to truncation error along the radial boundaries 
of the shearing box in our implementation of the shearing box boundary conditions, which 
remaps the EMFs rather than the magnetic field. We do this by gathering s = V ■ -B 



^ A Von-Ncumann Richtmyer artificial viscosity is a pressure proportional to (V • v)'^ that is small outside 
shocks. It should not be confused with the "a nomalous viscosity" used to model turbulent angular mom entum 
diffusion in standard accretion disk theories ( Shakura fc Sunvaevlll973t iLvnden-Bell fc Pringldll974[ ). 
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onto a single processor, solving the Poisson equation V^?/' = s using a standard FFT-based 
procedure, and then setting B B — 'Vip. The additional computational cost is negligible 
(usually less th a n 0.1 % of the total cost) because the operation is performed infrequently. 



Johansen et al.l (120081 ) have discussed other techniques to eliminate the rad ial dependence 



of num erical diffusion. This procedure eliminates the features reported in iJohnson et al. 

tosh . 



3. Turbulent Resistivity 

Our focus is on the diffusive effects of MHD turbulence induced by the magnetorotational 



instability (MRI; iBalbus fc Hawleylll99ll ). All models in this section have a mean toroidal 



field {B) = Boy, where Bo is chosen so that the initial plasma parameter f3 = SttPo/Bq = 
400. The models are evolved long enough (> 40 orbits) to reach a saturated (statistically 
steady) state. 

First we consider radial diffusion of a vertical field. Although this is the subdominant 
term in eq.([T]), it is easier to measure (for reasons we will discuss shortly) and thus allows 
us to explore systematic effects related to resolution, scale of the perturbation, and details 
of the initial state more readily. 

To measure the turbulent resistivity we evolve an initial state containing a mean field 
for hundreds of orbits. We then (arbitrarily) select an instant to inject a magnetic field 
perturbation of the form 

6Bz = asin{kj.x)z. (4) 

Here = 2T[nx/ is the radial wavenumber for the perturbation. We consider models with 
4 < Lx/ H < 32. We choose a so that it is larger than the background turbulent fluctuations 
Ifl, but small enough so that it does not greatly affect the background turbulence. Here we 
use 0.1 < a/Bo < 0.8. We then evolve the perturbed turbulence self-consistently, taking the 
sine transform of B^ to obtain a{t). 

A typical evolution a{t) is shown in Figure 1. The initial decay is approximately ex- 
ponential. We measure a decay rate 1/r by fitting an exponential to a{t) over at least one 
e-folding time. The decay rates are listed in Table 1. To give a sense of the amplitude of 
fluctuations in the background state we also plot in Figure 1 the evolution of the cosine 



•^The background contains power in the magnetic field a t small wavenumb er, with the power spectrum 
SB^ ~ const. ~ {SB'^)X^, where A is a correlation length (see Guan et al. 20081 for a discussion). The implies 
that the Fourier series coefficients in the background state will have rms amplitude {SB^)^^'^{X'^ / {LxLyLzY^'^ 
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amplitude; since the cosine amplitude remains so much smaller than the sine amplitude it is 

clear that the decay of a{t) is not simply due to phase drift. Eventually a{t) decays to the 
background level. In every case we have examined the perturbation decayed; the turbulence 
was stable to a large-scale magnetic field perturbation. 

Before going on we need to determine how strongly our measurement of r depends on our 
selection of initial conditions. We therefore measure r for an ensemble of initial conditions 
selected from the same run at widely separated times. We have done this for three models 
(s2, s5, and nl). The variation in r is < 13%, which may then be regarded as an error bar 
on our measurement. 

The decay time r is related to tjt as follows. Solving the induction equation with a 
scalar resistivity, 

a ^ aoexp{-r]Tklt), (5) 

so we define 

VT^{klr)-\ (6) 

Values of rj^ are given in Table 1. In general rjT will depend on the magnitude and direction 
of k and on the background field (B). 

Does the decay time scale as A;~^, i.e. does ifr depend on the magnitude of fe? In a model 
with (L-E, Ly, Lz) — (8, Att, 2)H we imposed perturbations with Ux = 1, 2, 3 on the same initial 
state. We find r = 32.7, 10.4,4.93 respectively, so that ?7t = 0.0495,0.0388,0.0366, crudely 
consistent with a diffusive scaling. Models nl, s5, and s6, with rij. = \ but — 8, 16, 32 and 
identical numerical resolution, have t/t = 0.0495, 0.0330, 0.0304, again crudely consistent 
with a diffusive scaling. So it appears that r^x is at most very weakly dependent on the 
magnitude of k. 

We checked the effect of resolution in models with N^l^x ranging from ?)2/ H to 128/ H. 
In all the runs Ax = Az ~ 2Ay (except in model s4, where Ax = Az = Ay; varying the 
zone aspect ratio made no difference in rj^) and {L^, Ly, Lz) = (4, 27r, 1)H. Comparing runs 
rl, r2, and r3, we find TjT = 0.037, 0.035 and 0.034 respectively. Our measurement of rjT 
is thus consistent with convergence, since the variation with resolution is smaller than the 
noise in decay time measurement that arise from choosing a particular initial state. 

We found that tit does depend on the field perturbation strength oq. In the {L^, Ly, Lz) = 
(8, 4:71, 2)H models with the same wavelength perturbation and a fixed resolution, when oq 
increases from O.ISq to QABq, rj^ increases shghtly (~ 20%) with oq; when oq = O.SBq tjt 
almost doubles. This is not surprising. It is known that the velocity fluctuation amplitude 
in shearing boxes increases in the presence of a background mean fleld. For sufficiently long 
wavelength sinusoidal perturbations, the imposed vertical flelds looks, locally, like a mean 
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field and so increases the velocity fluctuation amplitude and therefore rjT- 

We did not flnd any dependence of t^t on box size. Comparing runs with L^/H = 4 
and Ly/H = 27r,47r,87r we found r/T — 0.03 in every case. Comparing runs with L^/H = 4 
to 32, we again see that r^T lies in a narrow range around 0.03. Notice that for a vertical 
fleld perturbation with A^^ < 4H the decay time is less than Q~^. 

But what of the dominant term in eq([T]), the vertical diffusion of radial fleld? It is 
unclear how to measure this in the shearing box because the radial fleld is always buried 
"in the disk", where it is continually sheared into azimuthal fleld by the background shear 
(in a stratifled disk the radial fleld is present above and below the disk, where the Alfven 
speed is high and the plasma will corotate along fleld lines). So instead we measured the 
vertical diffusion of azimuthal fleld. This is also relevant to wind models because the radial 
fleld above and below the disk is always accompanied by an azimuthal fleld. 

To measure the azimuthal fleld diffusion we inject a perturbation of the form 

6By = asm{kzz)y (7) 
into an already turbulent state, measure the decay time r, and set 

VT = {klT)-\ (8) 

We imagine a large scale fleld entering the disk at a inclination of ~ 30 deg, running vertically 
through the midplane, and leaving at a similar inclination; to mimic this geometry we set 
2Ti/k, = L, = AH and 2H. 

The perturbation amplitude must be chosen larger than the turbulent background but 
as small as possible so it does not influence the background state. Because the turbulent 
fluctuations in the azimuthal fleld are larger than the fluctuations in the vertical fleld, a in 
eq.([7]) must be an order of magnitude larger (in comparison to Bq) than a in eq.(jl]). 

Model vl and v2 (see Table 1) have {L^, Ly, L^) = (4, An, A)H and resolution 128 x 200 x 
128; vl has ao = 2Bq while v2 has aoABo. We found 1/r = 0.052 and 0.050 respectively. 
This corresponds to a vertical diffusion coefficient //t ~ 0.020. Details are given in Table 1. 

To test for resolution dependence we repeated the above experiment at resolution 256 x 
400 X 256, with oq = ABq (model v3 in Table 1). We found 1/r = 0.040 and thus r/T ~ 0.016. 
This diffusion coefficient is ~ 20% smaller than that obtained at lower resolution. This 
small decrease in rj^ is mainly caused by slightly lower turbulent saturation level in the high 
resolution run in the perturbed state, where the saturation a is about 25% smaller. This 
might be surprising because on average higher resolution runs have higher saturation levels 
(see, e.g., Guan et al. 2009), but a fluctuates in time; our a is averaged over the same time 
interval used to flt 1/r. 
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4. Discussion: Turbulent Magnetic Prandtl number 

To calculate PrM,T, we need to assign a "viscosity" to the turbulence. We do this by 
measuring the turbulent shear stress 

Wxy,T = {pVxSvy - (9) 
and equating this to the shear stress that would be measured in a viscous fluid 

Wxy,v = pi^q^, (10) 
where q = —{l/2)dlnQ'^/dlnR = 3/2 for a Keplerian potential. Thus 

and the turbulent magnetic Prandtl number 

PrM,T-^ = (rl))(A:/ff^. (12) 
r]T q{p)ci 

The turbulent shear stress Wxy,T is related to a by a = Wxy^x/ {p)c^- The evolution of a from 
one of our runs is shown in Figure 1. 

We measure the time average of a during the decay time, denoted a (over the same 
period we fit r), and use equation ( |T2l) to calculate PrM,T- The results are compiled in Table 
1. The vertical field experiments give 0.35 < PrM,T < 0.58. The azimuthal field experiments 
give PrM,T ~ 1- 

In the current limited set of simulations we see no clear sign of PrM,T scaling with model 
or numerical parameters. This consistency is remarkable when we look at the radial diffusion 
of vertical field with different perturbation field strength qq. In the [L^, Ly, L^) = (8, An, 2)H 
models when the perturbation amplitude is strong, as in the ao = 0.8-Bo case, both r] and 
u!xy,T double compared to their weakly perturbed counterparts with ag = 0.1i?o; this doubling 
is precisely what is required for a constant PrM,T- 

We have also carried out comparison experiment with slightly larger initia l toroidal field 



strength Bo (0o = 100; model b2a in Table 1). Past numerical experiments (IHawley et al. 



19951 : iGuan et al.l l2008l ) imply that the saturation level scales linearly with Bq. For this 
model we found a = 0.0568 and r/x = 0.0878. Both the turbulent saturation level and r/x 
double, giving PrM,T ~ 0.43. 



Recently iLesur fc Longarettil (120081 ) have measured the turbulent resistivity in shearing 



box simulations. Their technique for measuring turbulent resistivity differs from ours: they 
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directly measure the EMF required to maintain a particular sinusoidal variation in the field, 
they measure a resistive stress tensor, and they consider only a mean vertical field (rather 
than the mean azimuthal field considered here). They find PrM,T = 1.6 for the diffusion of 
a radially varying azimuthal field. 

Transport properties of the MRI-generated turbule nt fiow have also been studied in the 
context of dust (passive sca l ar) mixing in a shearing box (Carballido et ahlboOS : Johansen fc Klahr 
20051 : 1 Johansen et al.ll2006l : iTurner et al.ll2006l : iFromang fc Papaloizoull2006l ). The analog of 
PrM,T in these experiments is the turbulent Schmidt number Sc = ut^/D^, where Dt is the 
diffusion coefficient for the grains. I n models with zero net magnetic flux (either when the 



dust is model ed as a passive scalar (ITurner et al. 



gas by drag ( Johansen &: Klahrl 



with a net vertical flux (jCarballido et al 



2005 



20061). or wh en the dust is coupled to the 
Fromang fc Papaloizoij 12006^) . Sc ~ 1. In models 



2005 



Johansen et al.l 120061 ) . Sc was found to be 



anisotropic and increase with a (and therefore depend on the initial vertical fleld strength), 
ranging from Sc ~ 1 for a weak fleld, to Sc ~ 2 for radial diffusion Sc ~ 10 for vertical 
diffusion when the mean fleld is strong and a ~ 0.5. Our results are broadly consistent 
with these measurements in the sense that Sc ~ PrM,T when the mean fleld is weak. Our 
scheme for measuring PrM,T is much more computationally expensive when the turbulence 
is strong, because a large computational volume is required to reduce the background fluc- 
tuations. It would be interesting to investigate whether PrM,T ^ 10 can be achieved with 
strong background flelds in future investigations. 



view. 



A PrM,T of ord er of unity is not surprising, perhaps, from a turbulent mixing point of 
Parked (119791 ). for example, argued that in isotropic turbulence i/t ~ ?7t ~ Iv, where / 



is the largest dimension of ed dies in the MHD tur bulent flow and v is the characteristic eddy 
turnover speed. Interestingly, lYousef et al.l (120031 ) also obtained an order of unity PrM,T from 
their turbulence simulations, where the turbulence is sustained by external forcing rather 
than the MRI in a disk. 

Can large scale magnetic flelds avoid escape from turbulent disks? We clearly have not 
included all of the effects outlined in the introduction that might influence the evolution of 
large scale flelds. Of these, perhaps the simplest route to large scale flelds is a rapid accretion 
mode in which the mean fleld torques the disk, causing it to accrete more rapidly than would 
a viscous disk. But our results cast doubt on models that conflne the large scale fleld by 
setting PrM,T ~ R/H (e.g. Shu et al. (2007)). 

Our models do not show a clear scaling of PrM,T with model parameters, but the 
dynamic range in parameters (and thus in ?7t) is small; future experiments over a broader 
range of initial fleld strengths may show scaling that is not evident here. Our models also 
do not include an explicit dissipation model, on which PrM,T might also depend. 
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Table 1. Model Parameters 
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Fig. 1. — Evolution of the amplitude a of a magnetic field perturbation SB^ — asm{kxx) 
in the presence of MHD turbulence (lower solid line; model s4). A linear fit to In a{t) is 
shown as a heavy dashed line. The decay time is r ~ 52J1~^, implying a turbulent resistivity 
r/T ~ 0.031. Also shown is the evolution of a in the same experiment (upper solid line). 
The turbulent magnetic Prandtl number PrM,T ~ 0.42. The dashed line shows the evolution 
of the corresponding cosine amplitude in SB^, indicating the amplitude of large scale field 
fluctuations in the turbulent background. 



